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Analysis of Extensive Cross-Flow Separation using Higher-Order RANS Closure Models 

J. H. Morrison * * * § A. G. Panaras t T. B.Gatski* G. A. Georgantopoulos § 


The turbulent flow fields associated with the incom- 
pressible flow over a 6: 1 prolate spheroid at high angle of 
attack, and the supersonic flow over an ogive cylinder are 
studied. Both these flows are characterized by large sep- 
aration and vortical flow regions and therefore provide a 
challenging database for comparison of turbulent closure 
models. Of interest is the ability to predict the effects 
of separation and associated vortical motion common to 
both flows. Two turbulent models are investigated that 
each represent the class of linear eddy-viscosity models 
(LEVMs) and explicit algebraic stress models (EASMs). 
Since the EASM accounts for anisotropic effects, the in- 
fluence of these effects on flow field predictions can be 
assessed. The EASM model is shown to both improve 
the separation location prediction and pressure trough 
under the secondary vortex on the 6:1 prolate spheroid 
at high angle of attack and high Reynolds number, and 
improve the prediction of the separation location on a 
supersonic ogive cylinder. 

Introduction 

Streamwise, or longitudinal, vortices are formed in 
various types of three-dimensional separated flows. The 
high-incidence flow about a slender body is the most 
widely known case. Crossflow separation occurs when 
fluid flowing circumferentially from the windward to the 
leeward side of such a body separates from the sides of 
the body along a separation line roughly parallel to its 
longitudinal axis. The fluid rolls up and forms two pri- 
mary vortices on the leeward side on both sides of the 
symmetry plane. The extent of the crossflow and the 
strength of the vortices grow as the angle of attack is 
increased. If the strength of the primary vortices is suffi- 
ciently large, secondary separation is induced below the 
primary ones, even if the boundary layer is turbulent. 
A similar type of separation is established in the high- 
incidence flow about delta wings, as well as in the flow 
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about various elements of high-speed vehicles (inlets, 
wing- and tail-fuselage junctures), where swept shock- 
waves interact with boundary layers and induce the ap- 
pearance of quasi-conical separation vortices. 

Though the numerical simulations of flows with exten- 
sive separation have significantly contributed to the un- 
derstanding of the nature of these flows, the accuracy of 
the predictions has not been exceptionally high. In var- 
ious Navier-Stokes calculations, that employ algebraic 
turbulence models and two-equation models, it has been 
found that the computations agree well with the data for 
moderate interaction strengths, but systematically under- 
predict the size of the interaction domain with increasing 
interaction strength. For example, in the case of high 
angle of attack flows, Panaras 1 reviewed calculations on 
a 6:1 prolate spheroid using algebraic and half-equation 
turbulence models. The review showed that algebraic 
models with appropriate modifications could predict sur- 
face pressures at high angles of attack relatively well; 
however, these algebraic models required user interven- 
tion and were not suitable for complex flows. This ad 
hoc modification of the algebraic turbulence model had 
the effect of reducing the eddy-viscosity coefficient in the 
region of crossflow separation. The same technique also 
proved successful in the case of swept shock turbulent 
boundary layer interactions. 2 A rationale for the success 
of these modified turbulence models that provide reduced 
turbulent stress levels within the separation vortices is 
easily provided. 3 When extensive cross-flow separation 
exists, the inner, turbulent part of the separated bound- 
ary layer winds around the core of the separation vortex. 
The outer layers, which are composed of low turbulence 
fluid, fold over the vortex and at the attachment region, 
penetrate into the separation bubble forming a low tur- 
bulence region that lies along the wall, underneath the 
vortex. The incorporation of this dynamic feature into 
a lower-order algebraic model, by replacing the equa- 
tion of calculation of the eddy-viscosity coefficient in the 
region of separation by a new one, which yields small 
values when the velocity gradients are high, brings con- 
siderable improvement to swept shock/boundary layer 
interactions, 4 crossing shocks and high angle of attack 
flows about slender bodies. 5 

Unfortunately, algebraic turbulence models suffer 
from their dependence on the definition of the length 
scale and are not viable for complex geometries. Higher- 
order turbulence models based on transport equations can 
account, in principle, for such transport effects as well 
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as turbulence anisotropies that algebraic models cannot. 
With this motivation, the performance of two-equation 
and explicit algebraic Reynolds stress models are eval- 
uated, using two high angle-of-attack flows with exten- 
sive cross-flow separation: an incompressible 6:1 prolate 
spheroid and a supersonic ogive cylinder. 

A variety of wind tunnel tests have been carried out 
for a 6: 1 prolate spheroid at DLR in Gottingen. The ma- 
jority of the tests were carried out at moderate Reynolds 
numbers, but some measurements were also conducted in 
collaboration with ONERA at higher Reynolds numbers. 
Ahn and Simpson 6 documented boundary layer transi- 
tion and separation characteristics for the 6:1 prolate 
spheroid for a range of Reynolds numbers and angles- 
of-attack. Chesnakas and Simpson 7 used an innovative, 
internally mounted LDV probe to measure the complete 
velocity vector and the Reynolds stress tensor from the 
viscous sublayer to the boundary layer edge at Re = 
4.2 x 10 6 . 

The review by Panaras 1 of calculations on a 6:1 pro- 
late spheroid were related to the DLR experiments. Tsai 
and Whitney 8 analyzed the prolate spheroid using a high 
Reynolds number K-e model with wall functions. Their 
results reproduced the flow structure and showed good 
agreement with measured pressure data at low Reynolds 
number conditions {Re = 4.2 x 10 6 ), but consistently 
under-predicted the measured windward pressure values. 
Kim et al. 9 investigated the prolate spheroid using two- 
equation and differential Reynolds stress models. For 
the low Reynolds number case (Re = 4.2 x 10 6 ) they 
matched the flow structure of the VPI results. However, 
the models that predicted the leeward pressures well typ- 
ically underpredicted the windward pressures, and the 
models that predicted the windward pressures well typi- 
cally overpredicted the leeward pressures. 

Constantinescu et al. 10 investigated the 6:1 prolate 
spheroid with both Reynolds averaged Navier-Stokes 
(RANS) and a detached eddy simulation (DES). DES is 
a hybrid approach that combines a RANS approach for 
near wall thin shear layers with a large eddy simulation 
(LES) of large scale turbulent structures. They pointed 
out that the boundary layer separation prediction is de- 
termined by the RANS model component of the DES 
procedure. They found that there were not substantial 
differences in the predicted values of wall pressure and 
skin friction between the RANS and DES approaches. 

An experimental study of the flow around a 3-caliber 
tangent ogive-cylinder in a supersonic flow was con- 
ducted at ONERA to develop an experimental database 
for code validations. 11 Panaras 5 investigated both these 
test cases with an algebraic turbulence model. The 
results at 30° angle-of-attack for the prolate spheroid 
match the experimental wall pressure distributions very 
well. In the case of the supersonic ogive cylinder, the 


agreement at 20° angle-of-attack was better than that at 
the lower angle (10°). These results indicated that addi- 
tional validation studies were warranted to further eval- 
uate the predictive performance of turbulence models in 
such high speed flows. 


Numerical Method and Turbulence Models 

The compressible Reynolds averaged equations for the 
conservation of mass and momentum can be written us- 
ing Favre-averaged variables as: 

+ (i) 


^ + (pUiUj+pdij) = - -j^r.PRi W 

where p is the density, Ui are the velocity components in 
the Xi coordinate directions, // is the molecular viscosity 
and patj is the mean stress tensor given by 


PtJ ij — p 


( dui duj\ 
\dxj dxi ) 


2 du k 
3dx k ij 


(3) 


The Reynolds stress tensor, , is a product of the av- 

eraging process and must be modeled to close the system 
of equations. 

Turbulence Models 

Chesnakas and Simpson 7 used a unique, internally 
mounted LDV system to measure the boundary layer on 
the 6:1 prolate spheroid. Their data showed that the flow 
gradient angle differs from the shear stress angle over 
most of the boundary layer. Linear eddy viscosity mod- 
els assume that the shear stress angle and the mean flow 
angle are the same. 

Rumsey and Gatski 12 investigated the effect of K- 
e and K-ui forms of two-equation models on multi- 
element airfoil. They demonstrated that models based on 
the K-e formulation were typically deficient in predict- 
ing adverse pressure gradient flows. Models based on the 
K -cu formulation provided a better prediction of adverse 
pressure gradient flows. Therefore, two models based on 
the K -oj formulation with and without non-linear terms 
are investigated in this study; the K-ui linear eddy vis- 
cosity model of Wilcox 13 and the explicit algebraic stress 
model of Rumsey and Gatski. 12 

K-ui Eddy Viscosity Model 

The K-ui model of Wilcox 13 is a linear eddy viscosity 
model that relates the Reynolds shear stress tensor to the 
mean strain rate tensor using the Boussinesq assumption: 


~PTij=Pt 




(4) 
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where the strain rate tensor is given by 
C _ 1 ( duj , 9tij \ 

u 2 \ dxj dxi ) 


(5) 


The eddy viscosity, // t , is related to mean and turbulent 
quantities by 

„ pK 

Pt = Cp — (6) 

w 

where K is the turbulent kinetic energy and ui is the spe- 
cific dissipation rate. 

The turbulent kinetic energy and specific dissipation 
rate are calculated by solving the following transport 
equations: 


^(pK) + -^(pUkK) =PV- C K2 pwK + 


d 

dx k 


P + 


Pt_\dK_ 


— (puj) + -^~{pU k U)) — Cui^V — C U 2PC0 2 + 


d 

dx k 

where the production is given as 

dui 


P + 


o K ) dx k 
— C u2 pU 
Pt \ 9 lj 


(7) 


o u ) dx k 


V = -r, 




OX; 


( 8 ) 


(9) 


and C, t = 1.0, Ck 2 = 0.09, C ul = 5/9, C u2 = 3/40, 
ok = 2.0, and <j w = 2.0. 


Explicit Algebraic Stress Model 

The explicit algebraic stress model of Rumsey and 
Gatski 12 replaces the linear Boussinesq approximation 
with the following non-linear relationship 


~PT%j — 2 /if { S tJ T [O'jtl-I ( Sj k 1 1 k j ^ 1 ' : ik S kj ) 

1 


- 2 ( 3 - 3 ^4 ( ' Sfcj k &ij 


H 


pK8 i:j (10) 


where 


Pt=~P Ka i (H) 

The strain rate tensor is given in Equation 5, the rotation 
rate tensor is given by 


Wa = 


1 

2 




( 12 ) 


and qi/t (t = 1/w) is obtained from the solution of the 
cubic equation: 

( t )' + '( t ) ,+ *( t ) + '“° <13 » 

See Rumsey and Gatski 12 for the definition of p , q, r and 
the solution of the above equation. 


The following K and ui transport equations are cou- 
pled with this non-linear stress-strain relationship: 


M (pK) + 


pV ~ fp*puK + 


d 

dx k 



dK' 

dx k _ 


(14) 


^ + 


d 

dx k 


f p* t \ duj ' 

[ fX+ ^)d^ k _ 


(15) 


where the coefficients are updated as in Rumsey 14 a k = 
1, (7<jj = k 2 / ~ 7)], k = 0.41, 7 = 0.53, 

B = 0.83, and C M = 0.0895. The function /,g* is from 
Wilcox 15 


where 


ft 3 * 


1 

l+680y| 

1+400 X ^. 


if Xk < 0 
if Xk > 0 


_ Cl dK di 0 
X" k w 3 dxj dxj 


(16) 


(17) 


Computational Approach 

The CFD code ISAAC 16 is used in this study. ISAAC 
is a second-order, upwind, finite-volume method where 
advection terms in the mean and turbulence equations are 
solved using Roe’s approximate Riemann solver 17 cou- 
pled with the MUSCL scheme. 18 Viscous terms are cal- 
culated with a central difference approximation. Mean 
and turbulence equations are solved coupled using an 
implicit spatially split, diagonalized approximate factor- 
ization solver. 

At the body surface, the no-slip boundary condition is 
imposed by setting the mean velocity component and the 
turbulent kinetic energy to zero, applying an adiabatic 
wall condition for the temperature, and setting the wall 
normal pressure gradient to zero. The wall condition for 
ui is due to Menter. 19 At inflow and farfield boundaries, 
conditions based on solving a one-dimensional Riemann 
problem are imposed on all solution variables. Symme- 
try conditions are applied at the symmetry planes and 
zero gradient conditions are applied at the outflow plane. 


Results and Discussion 

Two turbulent flow fields, one incompressible and the 
other compressible, characterized by large separation 
and vortical flow regions are investigated. The incom- 
pressible case is the 6:1 prolate spheroid at high angle- 
of-attack and high Reynolds number. The compress- 
ible case is the supersonic flow over an ogive cylinder. 
Both cases are predicted using the turbulence models 
described previously and results compared with experi- 
mental measurements. 
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Prolate Spheroid 

The prolate spheroid is a simple geometry that can pro- 
vide significant understanding of the complex flow that 
can develop about slender bodies at incidence. An at- 
tached boundary layer is formed on the windward side 
of the prolate spheroid. The flow accelerates around the 
body until it separates due to the adverse pressure gra- 
dient. The flow then rolls up into a separated vortex 
inducing a secondary flow separation and vortex roll up 
between the primary vortex and the body. 

As noted previously, a variety of wind tunnel tests 
have been carried out for a 6:1 prolate spheroid at DLR 
using a glass-fiber model of 4mm wall thickness having a 
major axis of 2.4m and a minor axis of 0.4m. Of interest 
here is the higher Reynolds number case tested at ON- 
ERA. The results presented here correspond to the 30° 
angle of attack. Re = 43 x 10 6 case. This case has pre- 
viously been used by an AGARD Working Group 20 to 
test boundary layer solvers. 

The experimental sting is modeled in the prolate 
spheroid calculations. On the surface of the spheroid, 
wall boundary conditions were imposed that specified 
zero velocity, adiabatic wall temperature, zero pressure 
gradient and zero turbulent kinetic energy. The calcula- 
tions were made on a half body using symmetry bound- 
ary conditions. At the outflow plane, all variables were 
specified with zero gradient, and characteristic based 
boundary conditions were used on the outer boundary. 
The grid extended approximately seven body lengths up- 
stream of the nose, six body lengths downstream of the 
end of the spheroid and approximately nine body lengths 
to the farfield. The experiment did not specify a transi- 
tion location. Transition is specified computationally at 
x/L k, 0.02, where L is the body length. 

A fine grid consisting of 129 points in the streamwise 
direction, 129 points in the wall normal direction and 73 
points in the circumferential direction (129 x 129 x 73) 
was generated using a two-dimensional hyperbolic grid 
generator and then rotating the grid about the axis of 
symmetry. A coarse grid (65 x 65 x 37) was generated 
by eliminating alternate grid points in each coordinate di- 
rection. The coarse grid had a maximum y + for the first 
point off the wall of approximately 0.67 with an average 
y + for the first point off the wall of approximately 0.44. 
A minimum of 25 points were located in the attached 
boundary layer regions. The fine grid had a maximum 
y + for the first point off the wall of approximately 0.3 1 
with an average y + for the first point off the wall of ap- 
proximately 0.21. A minimum of 60 points were located 
in the attached boundary layer regions. 

The solutions were computed at a Mach number of 
0.25. Other calculations were performed at nearby Mach 
numbers and showed no differences in the computed 
quantities. Solution residuals, surface pressures and skin 


friction were monitored to ensure iterative convergence. 

Figure 1 shows the top and side view of the surface 
streamlines calculated on the fine grid with the K-ui 
model. The topology of the flow is easily seen in this 
figure. The primary separation line forms at x/ L k, 0.1 
at the leeward symmetry plane and progresses windward 
down the body. The primary reattachment occurs in 
the plane of symmetry. Secondary separation and reat- 
tachment lines show the roll up of the secondary vortex 
between the primary vortex and the body. 



(a) Top View 



(b) Side View 

Fig. 1 Surface streamlines for prolate spheroid at 30° an- 
gle of attack calculated with the K-lo model on the fine grid. 

Surface pressure coefficients calculated on the coarse 
and fine grids for three streamwise locations are plotted 
as a function of circumferential location (<fi = 0 corre- 
sponds to the windward symmetry plane and 6 = 180 
corresponds to the leeward symmetry plane) in Figure 2 
for the K-ui model and in Figure 3 for the EASM model. 
The solutions on the coarse and fine grids for both the 
K-ui and EASM models at the x/L = 0.09 location in- 
dicate attached flow from the windward to the leeward 
side. This location is just upstream of the primary sepa- 
ration location shown in Figure 1 . The solutions for both 
the K-ui and EASM model show very little change at the 
x/L = 0.09 location with grid refinement. 

The solutions at x/L = 0.51 show a region of attached 
flow stretching from the windward plane to cf> m 130°. 
The attached flow region for the x/L = 0.81 location 
extends to cf> k, 110°. The coarse and fine grid solutions 
compare well with each other and the experimental data 
near the windward plane. However, as the separation lo- 
cation is approached, large variations arise in the coarse 
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and fine solutions for both the K-ui and EASM mod- 
els. The fine grid predicts an earlier separation than the 
coarse grid which is in better agreement with the experi- 
mental data. Additionally, the fine grid does a better job 
of predicting the pressure drop under the secondary vor- 
tex. Nevertheless, the variations in the pressure clearly 
indicate that the fine grid is still inadequate to resolve the 
flow. 




Fig. 2 Pressure coefficient for K-uj model on coarse and 
fine grid at x/L = 0.09, 0.51 and 0.81. 


The pressure coefficients calculated with the K -ui and 
the EASM model are shown in Figure 4 for the 129 x 





Fig. 3 Pressure coefficient for EASM model on coarse and 
fine grid at x/L = 0.09, 0.51 and 0.81. 


129 x 73 grid. Both models perform well at the upstream 
(x/L = 0.09) plane which contains only attached flow. 
The two models are virtually indistinguishable from each 
other and match experimental pressures very well. 

Both models predict the attached flow region at x/L = 
0.51 very well and agree well with each other. How- 
ever, the EASM model separates earlier than the K-ui 
model and predicts a stronger pressure trough under the 
secondary vortex. The models agree very well with 
each other at x/L = 0.81 only up to <f> « 40°; well 
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short of separation which occurs at <j> « 110° for the 
EASM model at this location. The pressures are con- 
sistently predicted lower than the experimental values at 
this streamwise location. The EASM model again pre- 
dicts an earlier separation and a larger pressure trough 
under the secondary vortex. The EASM model pre- 
dicts an earlier separation, a lower pressure plateau and 
a stronger pressure trough under the secondary vortex 
than the isotropic eddy viscosity model and more closely 
matches the experimental data. 



Fig. 4 Pressure coefficient for K-oj and EASM models at 
x/L = 0.09, 0.51 and 0.81 on fine grid. 


The primary separation line is shown in Figure 5 for 
both models on both grids. This figure confirms the re- 
sults shown in the pressure plots; the EASM model pre- 
dicts earlier separation on each grid than the K-ui model 
predicts. Figure 5 also confirms that the solution on the 
fine grid has not achieved grid independence and a finer 
grid is required. 



x/L 

Fig. 5 Separation line location predictions for K-ui and 
EASM models on fine and coarse grid. 

While the pressure predictions in the attached bound- 
ary layer regions of the prolate spheroid agree well be- 
tween the two models and the experimental data, the 
details of the predicted boundary layer differ substan- 
tially. Figure 6 shows a typical profile of the turbulent 
kinetic energy at x/L = 0.524 and <f> = 53.75° calcu- 
lated on the fine grid with each model (n is the distance 
normal to the surface). The EASM model clearly pre- 
dicts a lower turbulent kinetic energy out to the edge of 
the boundary layer than the K-ui model. To further illus- 
trate this point, the maximum turbulent kinetic energy in 
the boundary layer was calculated at each <f> location at 
the same streamwise location and is shown in Figure 7. 
The minimum in the turbulent kinetic energy maximum 
correlates well with the separation location. The EASM 
model consistently predicts a lower maximum turbulent 
kinetic energy through the attached boundary layer re- 
gion but then predicts a higher turbulence kinetic energy 
through the primary and secondary vortex. 

Although the results presented represent a self- 
consistent validation of two turbulence models for the 
flow over a prolate spheroid, it is apparent that the results 
still suffer from a lack of adequate resolution of the flow 
field. Further refinement in the grid structure needs to 
be undertaken in order to accurately predict the attached 
flow portion of this flowfield which is critical to accu- 
rately predicting the onset of separation. 

Supersonic Ogive Cylinder 

The supersonic flow about a 3 -caliber tangent ogive- 
cylinder at Mach 2.0 and angles of attack ranging from 
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Fig. 6 Turbulent kinetic energy profiles for K-uj and 
EASM models at x/L = 0.524, <f> = 53.75°. 



<t> 

Fig. 7 Maximum turbulent kinetic energy for K-lo and 
EASM models at x/L — 0.524. 

5° to 20° was investigated at ONERA. 11 The Reynolds 
number based on the body length was 1.44 x 10 6 . Exper- 
iments were performed for natural transition and transi- 
tion fixed at x/d = 1, where d is the body diameter. 

Grids were generated with a two-dimensional hyper- 
bolic grid generator by rotating the grid about the axis of 
symmetry. In total, three grids were utilized. A coarse 
grid consisting of 67 points in the streamwise direction, 
43 points in the wall normal direction and 37 points cir- 
cumferentially (67 x 43 x 37) and a medium grid con- 
sisting of 67 points in the streamwise direction, 85 points 
in the wall normal direction and 7 1 points circumferen- 
tially (67 x 85 x 71). Note that the streamwise spacing 
was maintained in order to accurately capture the bow 
shock. An additional fine grid consisting of 99 points in 
the streamwise directions, 115 points normal to the wall 
and 101 points circumferentially was used to investigate 
grid convergence on the a = 20° case. 

Two cases for the ogive cylinder were calculated; 
a = 10° and a = 20°. The 10° angle-of-attack case 
corresponds to moderately separated flow while the 20° 


angle-of-attack case is characterized by a large separated 
region. 

Inflow values of all flow quantities were held fixed at 
freestream values. Wall conditions were imposed with 
zero velocity, zero turbulent kinetic energy, adiabatic 
wall temperature and zero pressure gradient on the sur- 
face of the ogive cylinder. Outflow variables were set 
with a zero gradient condition. Symmetry conditions 
were imposed on the symmetry planes. Transition was 
specified at the experimental location, x/d= 1.0. 

The wall pressures predicted by the K-ui and EASM 
models on the 67 x 85 x 71 grid are shown in Figure 8 for 
the 10° angle-of-attack case at four streamwise locations. 
As in the prolate spheroid, <f> is measured from the wind- 
ward symmetry plane to the leeward symmetry plane. 
The agreement is good, especially at the downstream sta- 
tions, where the secondary separation is reproduced well 
by the calculations. The K-ui model consistently pre- 
dicts a lower pressure for the minimum pressure. The 
EASM model predicts a minimum pressure that is closer 
to the experimentally measured value. Additionally, the 
EASM model predicts a lower pressure under the pri- 
mary vortex (cf> « 155° at x/d = 7.5 and 8.5) that is in 
closer agreement to the experiment than the K-ui model. 

Figure 9 shows the grid convergence of the EASM 
model for the 20° angle-of-attack case. The separation 
location moves windward with increasing grid density 
and the pressure plateau in the separated region moves 
closer to the experimental value. The coarse grid is 
clearly insufficient to resolve this flow. However, there 
are small but still noticeable differences between the 
medium and fine grids especially in the separated re- 
gion and under the primary and secondary vortices. Even 
with the finest grid resolution the solution is not grid con- 
verged. 

The wall pressures predicted by the K-ui and EASM 
models on the 67 x 85 x 71 grid are shown in Figure 10 for 
the 20° angle-of-attack case. The agreement with data is 
very good except near (j) = 80° where the flow separates. 
The EASM model again predicts earlier separation and 
a pressure level at separation that is closer to the experi- 
mentally measured value. The K-ui model more closely 
matches the experimental pressures on the leeward side 
of the ogive cylinder at the two upstream locations. 

Concluding Remarks 

Two flows with extensive cross-flow separation were 
investigated with the Reynolds averaged Navier-Stokes 
equations using a two-equation eddy-viscosity model 
and an explicit algebraic Reynolds stress model. The 
flows were the incompressible 6; 1 prolate spheroid at 30° 
angle of attack and the supersonic (M = 2) ogive cylin- 
der at 10° and 20° angles of attack. It was found from 
grid convergence tests on both test cases that grid con- 
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Fig. 8 Pressure coefficient on ogive cylinder at 10° angle 
of attack for K-to and EASM models at x/d = 4.95, 6.0, 7.5 
and 8.5. 


vergence had not been achieved even with a 1 .2 million 
node grid for the prolate spheroid. Additional computa- 
tions on a refined grid are required to demonstrate grid 
convergence and to unambiguously evaluate the predic- 
tive capabilities of both models. 

The K-ui model of Wilcox 13 accurately predicted the 
surface pressures for the ogive cylinder at the lower angle 
of attack that had moderate separation. However, for the 
two cases with larger separations, the prolate spheroid at 
30° and the ogive cylinder at 20°, the K-ui model con- 






Fig. 9 Grid convergence of pressure coefficient on ogive 
cylinder at 20° angle of attack for EASM model at x/d = 
4.95, 6.0, 7.5 and 8.5. 


sistently predicted a smaller separation occurring later 
than the data showed. The EASM model of Rumsey and 
Gatski 12 consistently predicted an earlier separation with 
a larger vortical flow region than the K -cu linear eddy vis- 
cosity model, and is in closer agreement with the experi- 
mental data. Recall that one of the key advantages of the 
algebraic stress model is that it accounts for the effect of 
Reynolds stress anisotropy; whereas the isotropic eddy 
viscosity model does not. This causes the linear eddy 
viscosity model to be more dissipative than the EASM. 
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Fig. 10 Pressure coefficient on ogive cylinder at 20° angle 
of attack for K-uj and EASM models at x/d = 4.95, 6.0, 7.5 
and 8.5 on medium grid. 


The results obtained here, which generally show the im- 
proved predictive capability of the EASM, suggest that 
this better accounting of the turbulent stress anisotropies 
more accurately predicts the effect of the turbulence on 
the embedded vortical regions that characterize both the 
incompressible prolate spheroid and the supersonic ogive 
cylinder flows. Future studies into more accurately as- 
sessing and quantifying these anisotropic effects appear 
warranted. 
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